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We use infrared images obtained by the Thermal Emission imaging System (THEMIS) instrument on- 
board Mars Odyssey to retrieve the optical depth of dust and water ice aerosols over mote than 3.5 
martian years between February 2002 (MY 25, U ~ 330 - } and December 2008 (MY 29. L s ~ 

These data provide an important bridge between earlier TES observations and recent observation* from 
Mars Express and Mars Reconnaissance Orb iter. An improvement to our earlier retrieval (Smith, MD., 
Handheld, j.L, Christensen, P.R., Richardson, M.L 2003. j. Ceophys. Res. 108, doi:i0:1029/2003]E002114j 
to Include atmospheric temperature information from THEMIS Band 10 observations leads to much 
improved retrievals during the largest dust storms. The new retrievals show moderate dust storm activity 
during Mars Years 26 and 27, although details of the strength and timing of dust storms is different 
from year to year. A planet-endrdmg dust storm event was observed during Mars Year 2S near Southern 
Hemisphere Summer solstice. A belt of low-latitude water ke douds was observed during the aphelion 
season during each year, Mars Years 26 through 29, The optical depth of water ice douds is somewhat 
higher in the THEMIS retrievals at - 5:00 PM local time than in the TES retrievals at -2:00 PM, suggestive 
of possible local rime variation of clouds, 
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1. Introduction 

The continued successful operation of the Thermal Emission 
imaging System (THEMIS) instrument on-board the Mars Odyssey 
spacecraft has allowed for the long-term monitoring of the Mars 
atmosphere for more than 35 martian years (Christensen et a!„ 
2003). The atmosphere has been studied using both the visible and 
the infrared portions of the THEMIS instrument. Thermal infrared 
observations made by THEMIS in nine spectral bands enable the 
retrieval of the spatial and seasonal variation of dust and water ice 
column- Integra ted aerosol optical depth, as well as surface tem- 
perature and a vertically-averaged atmospheric temperature (Smith 
et aL, 2003), THEMIS visible images have been used to charac- 
terize the morphology and seasonal dependence of dust and iee 
douds (Jnada et ah, 2007), and to characterize mesoscaie douds 
(MtConnochie et ah, 2006). 

The retrievals of aerosol optical depth from THEMIS observa- 
tions provide a crucial link between the observations of aerosols 
taken by the Mars Global Surveyor TES (e.g. Smith, 2004) and Mars 
Orbiter Camera {MOO instruments (e.g. Cantor et aL, 2001: Cantor, 
2007) and those taken by instruments on-board the currently- 
operating Mars Reconnaissance Orbiter (MRG) (e.g. Wolff et aL, 
2009) and Mars Express (MEx) orbiter (e,g. Zasova et al, r 2005; 
Rannou et aL, 2006). The long overlap in time (more than one 
martian year) between the THEMIS observations and both MGS 
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and MRO/MEx observations provides an opportunity for the vali- 
dation and cross-calibration between all the different instruments. 
Together, these datasets now cover a period of more than five mar- 
tian years, enabling the initial exploration of interamiual variations 
in atmospheric conditions. 

The retrieval of dust and water ice aerosol optical depth using 
THEMIS infrared observations was previously described by Smith 
et a!. (2003), and more information on the THEMIS instrument 
can be found: in Christensen et a). (2003). A recent review of ah 
previous spacecraft observations of aerosols in the martian atmo- 
sphere is also given by Smith (2008). In this paper, we describe 
improvements to the Smith et aL (2003) retrieval as well as re- 
sults updated to the current time. In Section 2 we describe the 
set of THEMIS observations used in the retrieval of aerosol optical 
depth. In Section 3 we describe the retrieval algorithm, focusing on 
the improvements that have been made to the Smith e£ aL (2003) 
algorithm. In Section 4 we present the results of the retrieval along 
with a discussion of results, and we summarize our findings in 
Section 5. 

2. Dataset 

2.7. THEhUS instrument 

The THEMIS instrument obtains images of Mars using one of 
two separate focal planes. One focal plane contains ten spectral fil- 
ters covering the thermal infrared between 6,5 and 15 ptn, while 
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Fig, !. The distribution m soosan (1*5 and latitude of the THEMIS retrievals of aerosol optical depth. There is nearly-continuous coverage for more than 3.5 martian years. 
Small gaps in coverage are caused by solar conjunction and spacecraft anomalies, 

Che other contains five spectral fitters at visible wavelengths be* 25, L $ ^3304 to the end of December 2008 (MY 29, 4 = 183 s ). 

tween 0,45 and 0.85 pm (Christensen et aL, 2003). In this work Although the number of available images per day varies some- 

we use only the images taken in the thermal infrared. The first what with time, the coverage is nearly continuous over these 3.5 

two THEMIS spectral bands are identical and are centered at martian years with only short interruptions caused by solar con- 

6.78 pm. Bands 3 through 10 are centered at 7.93, 8.56, 9.35 t 10.21, junction or spacecraft anomalies. 

11.04, 11,79, 12,57 t and 14.88 pm, respectively with spectral widths Fig. 1 shows the seasonal (4) and latitudinal coverage of the 

that allow for slight overlaps between adjacent bands (except for THEMIS infrared images used in this work. Although the THEMIS 

Band 10, which is somewhat more narrow and has no spectral instrument cannot provide a continuous pole-to-pote mapping be* 

overlap with any other band). cause of data rate constraints, the coverage is still more than sttf- 

The THEMIS images are 320 pixels across-track and are built up ficient to provide a valuable picture of the seasonal and spatial 

in the along-track direction over time by spacecraft motion. Images trends of atmospheric aerosol optical depth, hi particular, the value 

are of variable length depending on the purpose of the observation of the atmospheric grids can easily be seen during the aphelion 

and the available resources available for downlink of the data. The season (4 = Q°~180 D } when data rates have typically been low. 
size of a THEMIS pixel at the surface is 100 m. The local time of the 

THEMIS observations varies between roughly 4:00 and 6:00 PM, 3, Revised retrieval algorithm 
which is somewhat later than either the 2:00 PM mean local time 

of Mars Global Surveyor or the 3:00 PM mean local time of die The retrieval used here is the same as Smith et af. (2003) algo- 
Mars Reconnaissance Orblter. rithm, except modified to include atmospheric temperature infor* 

We are most interested in the variations of atmospheric con* mation from the THEMIS observations themselves. Here we give a 

ditions on scales larger than, the THEMIS pixel size, so for each brief outline of Smith et al, (2003) algorithm, and in the next sec- 

THEMIS image we use as our observation the average of a block of bon we describe the revision to include temperature information 

data 320 pixels wide (the entire width of the image) by 256 pixels from the THEMIS observations. 

long. The last 256 pixels of each image are used for this average The basic idea of the THEMIS retrieval is to find the values for 
because they have the most accurate calibration (an effect only no* surface temperature and dust and water ice aerosol optical depth 

ticeable for Band 10). The resulting average covers an area 32 by that provide the best fit between the computed and observed ra- 

26 km in size, or roughly one*half degree square, diance spectra in THEMiS Bands 3-8 (roughly 8-12 pm). The three 

quantities, surface temperature, dust optical depth, and water ice 
2.2. Observations used in this study optical depth are fit simultaneously by linearizing rhe solution for 

radiance about the current best guess and iterating until solution. 

THEMIS images are divided between those targeted at objects Convergence typically occurs in less than five iterations, 
of geologic interest, general mapping of the entire surface, and a Radiance is computed using a purely absorbing plane-parallel 

latitude-longitude grid of observations taken every couple weeks atmosphere. Thus, the optical depth values reported here should 

intended to survey atmospheric conditions. Although the targeted be considered as effective absorptive optical depth only, not full 

observations generally provide a reasonable global-scale sampling extinction optical depth that includes scattering, Numerical exper- 

over time, the atmospheric grid has proven to be very useful as iments show that extinction optical depth is roughly 1,3 times 

a means to "fill iirt any spatial gaps in the coverage of targeted absorption optica! depth for dust, and 1.5 times absorption opri- 

observations. cal depth for water Ice aerosol (Smith, 2004). 

Included in this work are all THEMIS images that include ail The spectra! dependence of optical depth for dust and water ice 

ten infrared bands and for which the surface is sufficiently warm aerosols is taken from TES retrievals (Bandfield and Smith, 2003) 

(220 K or higher) to provide enough thermal contrast between and are assumed to be constant in space and time. Because surface 

surface and atmosphere to allow for a reliable retrieval. A total emissivity has a spectral dependence that is very similar to that of 

of 26,400 images fit these criteria, covering the period from the dust aerosol over the THEMIS wavelengths, the two cannot be re- 
beginning of the mission in February 2002, or Mars Year (MY) frieved independently. Instead, TES results are used for the spectral 
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Normalized Contribution Function 



Fig. 2. Two detertninsUoifcs of the equivalent height - represented by THEMIS Band 10 temperatures. (Left) Computed normalized contribution functions for a temperature profile 
typical of daytime conditions at low latitudes and U = 180° t and for two isothermal atmospheres bracketing expected martian conditions. (Right) Observed temperature 
differences between concurrent THEMIS Band 10 and TES retrieved temperatures at different pressure levels. Both methods indicate that THEMIS Band 1G is representative 
of a wide range of heights centered at about G,$ mhar. 


dependence of surface emissivity (Handheld and Smith, 2003) and 
for the amplitude of surface emissMty as a function of latitude 
and longitude {Smith, 2004). 

The vertical distribution of dust aerosol is assumed to be well- 
mixed. Water ice aerosol is assumed to be in the form of conden- 
sate clouds, so that there is no aerosol below the water condensa- 
tion level The condensation level is computed using a water vapor 
abundance taken from concurrent TES observations, or from TES 
climatology after the end of systematic TES spectrometer observa- 
lions (MY 27, 

3 A. Inclusion of THEMIS temperature information 

The one THEMIS band sensitive to atmospheric temperature 
(Band TO) was designed to give a temperature representative of 
a broad range of the atmosphere (roughly 0.1 ~2,0 mbar} ¥ simi- 
lar to the *T15 H temperature derived from Viking Orbiter IRTM 
observations (Martin and Kieffer, 1979). Although useful this sin- 
gle vertically-integrated temperature by itself is not sufficient to 
accurately model atmospheric radiance. The Smith et al (2003) al- 
gorithm relied on concurrent retrievals from TES to provide the 
necessary temperature profiles. That strategy is no longer possible 
after the failure of the TES spectrometer (near-continuous obser- 
vations by the TES spectrometer ended on 31 August 2004, or 
MY 27 f I* =■ 81*). For more recent THEMIS observations we use 
a combination of the averaged atmospheric temperature provided 
by THEMIS Band 10 with historic observations by TES at the same 
season and location (i.e. climatology) to estimate the temperature 
profile, T(p), al the time of each THEMIS observation. 

Although the temperature derived from THEMIS Band 10 is rep- 
resentative of a broad vertical range of the atmosphere, it is useful 
to define a single effective height for this temperature. Fig. 2 shows 
both a theoretical and an empirical derivation of this height. The 
left panel of Fig. 2 shows the contribution function for THEMIS 
Band 10 computed by taking the vertical derivative of the trans- 
mittance of the CO 2 atmosphere as a function of height integrated 
over the spectral response of THEMIS Band 10, The contribution 
function shows the vertical distribution of the relative contribu- 


tion to the observed radiance in THEMIS Band 10. This distribution 
is not sensitive to the temperature profile as shown by the very 
similar contribution functions for three widely different tempera- 
ture profiles. The peak of the contribution function is at a pressure 
level of about 0.6 mbar, although the peak is broad. The right panel 
of Fig. 2 shows the rms temperature difference between observed 
THEMIS Band 10 brightness temperatures and the temperatures at 
different pressure levels extracted from concurrent TES retrievals. 
This comparison includes ail THEMIS observations taken during the 
period when the TES spectrometer was still operational, and so in- 
cludes a complete range of seasons and latitudes. The minimum 
temperature difference is found at a pressure level of 0.5 mbar, 
which is in reasonable agreement with the theoretical results given 
the wide vertical distribution apparent from both methods. For the 
purpose of this retrieval we assign an effective height of 0.5 mbar 
to the THEMIS retrievals. 

Much of the time, using TES climatology temperatures for the 
THEMIS retrieval is a reasonable approximation, Fig. 3 shows a 
comparison between TES climatology temperatures at 0,5 mbar 
and observed THEMIS Band 10 temperatures. The TES climatol- 
ogy is taken from the first martian year of TES operations, from 
MY 24, U = 104* to MY 25, U =.• 104 s5 . There is excellent corre- 
spondence between the two during the aphelion season, L s ™ 0° 
to ISO*. However, it is also apparent that the random timing and 
strength of large regional and planetary-scale dust storms during 
the perihelion season, L — 180° to 360 c , makes the use of TES cli- 
matology at those times much less desirable. These temperature 
differences can be quite large, exceeding 30 K for the largest dust 
storms. When actual atmospheric temperatures are greater than 
those predicted by TES climatology, then the thermal contrast be- 
tween the surface and atmosphere will be overestimated, which 
means that aerosol optical depth will be underestimated. 

It was noted by Smith (2004) that the atmospheric warming 
caused by dust storms varies with height. The largest response is 
at high altitudes, while there is little response at levels near the 
surface. Fig, 4 shows a time series of temperatures at low Southern 
latitudes from TES over three martian years for five different atmo- 
spheric levels. While there is large interannual variation from year 
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Fig. 3. A comparison of observed THEMIS Band 10 temperature against those estimated using TES climatology (at 0.3 mbar} from the first year of TES observations (MV 24 
and 25}. TES climatology provides a good estimate except when there are large dust storms. 



Fig. 4, Atmospheric temperatures retrieved from TES data as a function of season (L=) and pressure level at a given location (20” 5. 270“ W). Variation m temperature from 
year to year in response to dust storms is most pronounced at greater heights above the surface and becomes relatively small near the surface 



Latitude 


Fig. S. The difference in zonaily-averaged thermal structure of the martian atmo- 
sphere shown as a cross-section in latitude and log-pressure (the 0;s~mbar pressure 
level is ^40 km above the 6.1 mbar level) as observed by ITS. The gently slop- 
rug line across the bottom indicates the zonaliy-averaged pressure at the surface. 
Shown is the change in the daytime temperatures as a result of the 2001 (MY 25) 
planet -encircling dust storm event This Is temperature at MY 25 T L* — 2G5 ,£> (July 
2001} minus the temperature at MY 24, L A -205“ (September 1999). The temper- 
ature difference rises almost linearly from near-zero at the surface to a maximum 
value somewhere above the O.S-mbar level. 

to year caused by dust storms at 05 mbar, the year-to-year vari- 
ation near the surface (e,g. ¥ 37 mbar) is much reduced, Another 
example is given by Fig. 5, which shows the temperature differ- 


ence as a function of latitude and height between temperatures 
observed by TES at MY 25, U = 205° at the peak of a planet- 
encircling dust storm and those observed exactly one martian year 
previously when there was no dust storm (Smith et al. + 2002). Here 
we see a near-zero temperature change at the surface increasing 
roughly linearly with height. 

The TES observations of temperatures during dust storms sug- 
gest a relatively simple correction based on THEMIS Band 10 tem- 
peratures that can be made to improve the estimate of the tem- 
perature profile from TES climatology. We define a quantity, AT, 
as the difference between the observed THEMIS Band 10 tempera- 
ture and the TES climatology temperature at 05 mbar (at the sea- 
son, latitude, and longitude of the THEMIS observation). Recalling 
Fig. 3, most of the rime AT is within 3 K of zero, but during dust 
storms AT becomes large and positive with values of 10-40 K, The 
temperature profile at the time of the THEMIS observation. Tip), 
is then estimated from the TES climatology temperature profile, 
Ttes(P). rising; 

T(p) = Ttes (P ) 4- 0.4 A r log — , ( 1 } 

P 

where p re f is a reference pressure level taken to be 6.1 mbar. Equa- 
tion (1) provides a temperature correction that is zero at the ref- 
erence pressure (essentially the surface) and increases linearly in 
log-pressure (nearly linear in height) to a value of AT at a pres- 
sure level of 05 mbar. 
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L s 


Pig, 6. An overview of THEMiS rerrieved aerosol optical depth and Band 10 temperature over more than 3.5 martian years, Shown is the zona! average of each quantity as a 
function of season (L s ) and latitude. (Top) Dust optical depth at 1075 cm - ' scaled to an equivalent g.l-mbar pressure surface to remove the effect of topography, (Middle) 
THEMIS Band 10 temperature, representative ef a wide range of heights centered at about 0.5 mbar. (Bottom) Water ice optical depth at S25 cm - ’ 1 , 


32. Estfmafion of uncertainties 

A number of sources contribute to uncertainties in the retrieved 
values of dust and water ice optical depth, including instrument 
noise and calibration, error in the assumed temperature profile and 
surface emissivity, and simplifying assumptions of fixed spectral 
shapes and a well-mixed non-scattering aerosol. The formal prop- 
agation of random instrument noise leads to negligible error since 
we average over 320 x 256 pixel frames. Uncertainty from calibra- 
tion and other systematic errors is not reduced by averaging pixels. 
Such uncertainties are difficult to evaluate, but our experience with 
the TES retrievals, and the comparison of THEMIS optical depth re- 
sults with TES leads to an estimate of an uncertainty for a single 
THEMIS observation of 0.04 or 10% of the total absorption opti- 
cal depth, whichever is larger, Uncertainties are likely somewhat 
higher (perhaps 20% or even greater) during the most intense dust 
storms because large corrections to the temperature profile must 
be made for those observations. The estimate for the uncertainty 
in optical depth also does not include the systematic difference 
between absorption and extinction optical depth. As stated earlier, 
we estimate that the extinction optical depth (including scattering) 
is systematically higher than the absorption optical depth reported 
here by ^30% for dust and —50% for water ice, 

4. Results 

The retrieval of dust and water ice aerosol optical depth was 
performed on all THEMIS daytime infrared images with surface 
temperature greater than 220 K and all ten spectral bands present. 


The images span the entire 3,5 martian years of Mars Odyssey op- 
erations from MY 25, L> ~ 33Q C (IS February 2002) to the present, 
MY 29, L s ~ 183* (31 December 2008 }, yielding a total of 26,400 
good retrievals, 

Fig, 6 shows the results of our retrievals. Because dust has 
been assumed to be well-mixed, in Fig. 6 and the following 
figures we present dust optical depth scaled to an equivalent 
6,1-mbar pressure surface to remove the effect of topography. This 
is accomplished by multiplying dust optical depth by the factor 
(6,t mbar/surface pressure). Water ice aerosols have not been as- 
sumed to be well-mixed, and so we present water ice optical 
depth without any scaling. Band 10 temperatures are shown for all 
daytime infrared images without the restriction on surface tem- 
perature. 

4.1. Dust optical depth 

The history of dust optical depth shows the familiar seasonal 
dependence observed by previous spacecraft (Smith, 2008). The 
highest dust optical depth on a global scale is found during the 
perihelion season (L s ~ lS0 4< -360°), but with significant variation 
in the timing and amplitude of large dust storms from one mar- 
tian year to the next Relatively low dust optical depth persists 
throughout most of the aphelion season with the exception of oc- 
casional small-scale dust storms at high latitudes along the edge 
of the retreating polar ice caps. 

The MY 26 dust storm season was relatively mild. As usual, 
dust activity began to pick up after l s = 140° with the first large 
regional storms at around 1* = 210°. Additional dust activity was 
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observed at high southern latitudes near summer solstice. The dust 
storms at Is ™ 315 s are a regular feature of the dust cycle, but were 
noticeably more intense than those seen in previous years at that 
season (Liu et al., 2003; Smith, 2004). These dust storms were par- 
ticularly notable since they occurred at the time of the Beagle 2 
arrival at Mars, and just prior to the landings of two Mars Explo- 
ration Rovers (MER). 

The MY 27 dust storm season was also mild, although there 
were several differences in the timing and intensity of the largest 
dust storm events compared with MY 26, An unusually early 
and intense low-latitude dust storm was observed in MY 27 at 
L s — 135°, which was also observed a£ both MER landing sites by 
the rovers (Smith et at, 2006). However, this early activity was 
not followed up by increased dust optical depth during the rest of 
the season. Regional dust storms at 4 = 225° were comparable in 
intensity although later than those observed in MY 26. The late- 
season dust activity at l s ~ 310° was more similar to the moderate 
dust storms observed in MY 24 and 25 by TES (Liu et al. 2003; 
Smith, 2004) than the more intense dust event of MY 26. 

By far, the greatest dust activity observed so far by THEMIS 
occurred during MY 28. The early-season dust activity was not ex- 
ceptional, with the regional dust activity at 1$ — 220° somewhat 
less intense than in previous years. However, at l s = 265° a se- 
ries of dust storms began, which lead to a planet-encircling dust 
event. Thermal infrared dust optical depth exceeded unity over 
most of the planet for several weeks making this the highest global 
dust loading since the MY 25 (2001) event (Smith et aL, 2002; 
Cantor, 2007), Dust optical depth only gradually decayed, staying 
above 0.5 until at least l s ^ 310 s . No new significant dust ac- 
tivity was observed by THEMIS during the usual l s ~ 3 10° -320° 
dust storm period, presumably because of the large perturbations 
of conditions from the late-season planet-encircling dust event. 

Observations from THEMIS as we enter the dusty season for MY 
29 once again show elevated eariy-season dust activity similar to 
that observed in MY 27. 

42. WEMtS Band 10 temperatures 

The middle pane! of Fig. 6 shows THEMIS Band 10 bright- 
ness temperature, which was shown above in Section 3.1 to be a 
broad vertical average of atmospheric temperature with an effec- 
tive height of roughly 0,5 mbar or 25 km above the surface. At- 
mospheric temperatures respond to a combination of orbital (per- 
ihelion vs. aphelion) and seasonal (summer vs. winter) variations, 
as well as to direct heating by the absorption of sunlight by dust 
and to genera! circulation patterns. The observed heating caused 
by the large regional storms during MY 26 and 27 was roughly 
15-20 K, while that caused by the planet-encircling dust event of 
MY 28 was 30-40 K. Globally-averaged temperatures have a con- 
sistent and well-defined minimum near l s — 40 s , which is earlier 
than both aphelion (L $ =71°) and the global minimum dust opti- 
cal depth (U ™ 130'), a phenomenon also observed by TES (Smith, 
2004). 

4.3. Water ice optical depth 

The main feature evident in the THEMIS retrievals of water ice 
aerosol optical depth (bottom panel of Fig. 6) is the aphelion sea- 
son low-latitude dcud belt. This robust set of cloud features has 
been observed repeatedly by telescopic observations fe.g. Clancy et 
ai. + 1996), Viking (e.g. Tamppari et aL, 2000), TES (Pearl et aL, 2001 ; 
Liu et al, 2003; Smith, 2004), and other instruments. Low-latitude 
water ice optical depth increases rapidly after L s — 0° reaching its 
peak value by about L = 6Q C . The cloud-belt begins to dissipate 
after L s = 140 s as atmospheric temperatures rise, but remnants 
persist until at least l s = 180°. At its peal<, there is significant cloud 



Fig. 7, A comparison of retrievals of dost optical depth from THEMIS observations 
using TBS climatology and using THEMIS Bond 10 information. The red points 
are those taken during the 2007 planet-encircling dust storm event (MY 28, t< — 
265" -3054. 

cover between 10 s S latitude and 30 c ' N latitude, with enhance- 
ments over areas of elevated topography such as Elysium, Tharsis, 
Olympus Mons, and Alba Patera. There is relatively little interan- 
rtuaf variation in the timing, location, and amplitude of the cloud 
belt However, the THEMIS observations do show that very early 
low-latitude dust storms, such as at MY 27, L s = 135° , can very 
quickly reduce cloud optical depth to near-zero values. In this par- 
ticular case, the cloud belt partially reformed by l s = )50 c after 
the dissipation of the dust storm and atmospheric conditions had 
returned to nominal values for that season. 

The polar hood clouds that form at high latitudes in the winter 
hemisphere are not observed by THEMIS because at the relatively 
late local time of the Mars Odyssey orbit (—5:00 PM) surface tem- 
peratures are too low at those latitudes to provide the thermal 
contrast between the surface and atmosphere necessary for the re- 
liable retrieval of water ice cloud optical depth. 

4.4. Comparison against retrieval using TES climatology 

The main difference between the original retrieval algorithm of 
aerosol optical depth from THEMIS infrared images used by Smith 
et al, (2003) and the algorithm used in this work is in how the 
atmospheric temperature profile is obtained. In the Smith et at. 
[2003) algorithm, concurrent TES observations were used to es- 
timate atmospheric temperatures, but that is no longer possible 
since the end of systematic TES spectrometer observations on 31 
August 2004 (MY 27, l s - 81°). Instead, we estimate atmospheric 
temperatures using a combination of observed THEMIS Band 10 
temperature and historic TES observations (climatology) from the 
season and location of interest. 

Fig. 7 shows dust optical depth retrieved using observed 
THEMIS Band 10 temperature information compared to that re- 
trieved using TES climatology alone without correction. In the 
figure, the red points are those observations taken during the 2007 
planet-encircling dust storm event (MY 28 T l s = 265 W -305 G ), while 
the black points are those taken at other times. The blue lines 
show when the two dust optical depth results are the same within 
uncertainties. 

As expected, when there are no large dust storms the difference 
between retrieved dust optical depth using TES climatology and 
the correction using THEMIS Band 10 temperatures is essentially 
within the uncertainty of the algorithm. However, during large 
dust storms it is necessary to use observed THEMIS Band 10 tem- 
peratures to get an accurate retrieval During large dust storms the 
observed THEMIS Band 10 temperature is higher than the climatol- 
ogy temperatures from the first year of TES observations (MY 24, 
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Ftg, 8. A comparison of dust anti water ice optical depth retrieved from THEMIS and TES f>bservatiops during the time when both instruments were operational t MV 2$, 
L ~ 330 s io M¥ 27, L s - 8F). The scale for dust optical depth is from 0 (purple) to 03 (red), The scale for water ice optical depth is from 0 (purple) to G.IS (red). 


l s ~ 104° to MY 25, 4 = 104°), which was a very mild year for 
dust storms (Smith, 2004}. Atmospheric temperatures greater than 
climatology lead to a smaller thermal contrast between the surface 
and atmosphere than vvould be expected from climatology. There- 
fore, to produce a spectral feature of a given size, a greater amount 
of dust must be present than would be required from atmospheric 
temperatures from climatology. 

Fig. 7 shows that the above effect can be quite significant dur- 
ing the largest dust storms: The dust optica! depth can be underes- 
timated by up to a factor of two if observed temperatures are not 
used. Dust optical depth can also be underestimated by a smaller 
amount (10-20%) during large regional storms as shown by the up- 
ward curve of black points at larger values of dust optical depth. 

4.5. Comparison against concurrent TBS observations 

Regular observations by the TES spectrometer continued for 
more than a martian year after the beginning of THEMIS obser- 
vations allowing for a direct comparison between the two instru- 
ments and a validation of results from THEMIS. Fig. 8 shows side- 


by-side plots of dust and water ice optical depth retrieved from 
TES and THEMIS. The correspondence between the two is close, 
with both the seasonal and spatial patterns and the amplitudes of 
results matching welt Both TES and THEMIS observe the decay of 
a regional dust storm at the end of MY 25, the usual low dust 
optical depth during the aphelion season, and a relatively mild 
dusty season in MY 26 with moderate regional storms at 4 = 210° 
and 315°. The aphelion season water ice cloud belt is observed by 
both instruments to begin intensification after l s = 0°, reach full 
strength at L s = 60° with the same amplitude and latitude extent, 
and linger until about L - 180 ! *. 

One difference is the apparent cloudier conditions observed 
by THEMIS between U ™ 120° and 180°. This was interpreted 
by Smith et ai. (20Q3) as a result of the later local time of the 
THEMIS observations. Only one year (MY 26) in that seasonal pe- 
riod was observed by both TES and THEMIS concurrently. How- 
ever, observations during this season for other years (MY 24 and 
25 for TES; MY 27, 28, and 29 for THEMIS) show that during 
roughly L s = tGG ,! ~18G° the water ice optical depth observed by 
THEMIS is always somewhat higher than that observed by TES, ex- 
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cept for the sudden drop in clouds observed by THEMIS at MY 27, 
l s - 135* at the onset of an unusually early low-latitude dust 
storm. The difference in optical depth, is about 0,03, which is dose 
to the level of uncertainty. Nevertheless, the consistency of the 
trend over several years of observations from both instruments 
gives support to the idea that cloud optical depth varies with local 
time, a tendency also observed from telescopes (Wolff et aL 1999; 
Akabane et aL 2002; Glenar et aL 2003), Viking Orbiter (Tamppari 
et aL, 2003 }, Mars Global Surveyor (Wilson et aL, 2007), and Mars 
Reconnaissance Orbiter (Maim et aL, 2008). 

Notable in Fig, 8 is the smaller latitudinal range given by the 
THEM15 observations because of its later local time and thus cooler 
surface temperatures. This prevents the THEMIS observations from 
observing the polar hood water ice clouds seen by TES as well 
as some of the smaller dust storms along the edge of the polar 
caps. At the present time the Mars Odyssey spacecraft is per- 
forming maneuvers to move the local time of the orbit earlier, 
to about 3:00 PM, The move is expected to be complete in late 
20O9. This should: provide a somewhat greater latitude extent for 
future aerosol optical depth retrievals, and may help provide a 
more definitive conclusion about the diurnal variation of water ice 
clouds in THEMIS observations. 

5. Summary 

Infrared images taken by the THEMIS instrument on-board Mars 
Odyssey are welt-suited for the retrieval of dust and water ice opti- 
cal depth in the Mars atmosphere. Use of atmospheric temperature 
information contained in THEMIS Band 10 data in the retrieval 
algorithm provides a much more accurate retrieval during large 
dust storms as compared to using temperatures from a climatology 
database alone. During the planet-encircling dust storm event of 
2007 (MY 28), the observed THEMIS Band 10 atmospheric temper- 
atures were 30 K warmer than those from TES climatology (Smith, 
2004). Using a temperature profile based on the observed THEMIS 
Band 10 temperatures instead of directly from TES climatology re- 
sulted in an increase in retrieved dust optical depth during the 
storm by a factor of two or more. 

The long and nearly-continuous record of THEMIS retrievals 
of aerosol optical depth provides a critical link between historic 
observations from Mars Global Surveyor data and the recent ob- 
servations by the Mars Express and Mars Reconnaissance Orbiter 
spacecraft. Comparisons of aerosol optical depth retrieved from 
THEMIS and concurrent TES observations shows strong agreement 
providing a validation of the consistency of results from the two 
instruments. Similar cross-correlation and validation is currently 
underway comparing THEMIS retrievals against those from Mars 
Express and Mars Reconnaissance Orbiter instruments (e.g. Smith 
et aL 2009; Wolff et ah. 2009). 

Over a time span of more than 3.5 martian years, THEMIS ob- 
servations have shown both repeatable patterns and large inter- 
armual variation in aerosol behavior, Mars Years 26 and 27 were 
observed to have moderate dust activity during the perihelion 
(southern spring and summer) season, with regional dust storms 
near L s =^220* and 315°, while MY 28 featured a planet-encircling 
dust storm event just before southern summer solstice but no late- 
season regional storm at L = 315'-. Early season low-latitude dust 
storms were observed in both MY 27 and 29 before L = 180 s , a 
type of event not observed by TES during MY 24-26. 

The low-latitude aphelion season water ice cloud belt was ob- 
served each year from MY 26-29 to begin forming at L s = 0T 
reaching peak optical depth between L s = 100* and 120°, and 
gradually decay after L s = 140°. Cloud thickness was observed to 
rapidly decline during the early-season dust storms in MY 27 and 
29, and to recover to nominal values after dissipation of the dust 
storms. 


The limited latitudinal extent of reliable THEMIS retrievals be- 
cause of the late local time (-5:00 PM) of the Mars Odyssey orbit 
precluded the observation of both polar hood clouds and small 
dust storms along the edge of the polar caps. However, the later 
local time of THEMIS observations compared to TES (-2:00 PM) 
does allow for a limited examination of the variation of cloud op- 
tical depth with time of day. The THEMIS observations consistently 
show somewhat higher water ice optical depth in the low-latitude 
aphelion season doud belt, especially after Is — 100*. 

The THEMIS instrument remains fully operational and continues 
to add to its long history of atmospheric observations. The Mars 
Odyssey spacecraft is currently undergoing maneuvers to bring the 
local time of its orbit to an earlier time. It is expected that the ma- 
neuvers will be completed in late 2009 with the orbit near 3:00 
PM local time, Observations at the earlier time will allow the re- 
trieval of aerosol optical depth over a greater latitude range and 
will provide further information on the ti me-of-day variation of 
water ice optical depth. 
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